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A SOLUTION OF THE NAVIER-STOKES EQUATION FOR SOURCE AND 
SINK FLOWS OF A VISCOUS HEAT-CONDUCTING 
COMPRESSIBLE FLUID 
By Robert V. Hess 

SUMMARY 


A solution of the Navier-Stokes equation for source and sink flows of 
a viscous, heat-conducting, compressible fluid is given for the case of 
constant total flow energy. For the satisfaction of the condition of 
constant total flow energy, a certain Prandtl number is required. Aside 
from this more obvious requirement, the selection of a certain ratio of 
the first and second viscosity coefficient is necessary. The nature of 
the general solutions for flow with arbitrary Prandtl number and with 
heat addition is discussed. Furthermore, the manner is discussed in which 
the familiar heat-conduction effects combined with the peculiar viscous 
effects solely due to compressibility, sometimes called the longitudinal 
viscous effects, influence the flow through a curved minimum section 
joined to a sink flow. A discussion of the second viscosity coefficient 
from the gas-dynamic approach is also given. 


INTRODUCTION 


With the advent of flight at extreme altitudes, the study of high- 
speed viscous flows has become of practical importance. Since for high- 
speed flows the compressibility effects are large, the viscous effects 
solely due to compressibility also called the longitudinal viscous effects 
have to be investigated in addition to the transverse viscous effects 
which are associated with boundary layers. 

Since the transverse viscous effects and the longitudinal viscous 
effects in one dimension are well-known, there remain to be investigated 
the longitudinal viscous effects in more than one dimension and their 
interaction with the transverse viscous effects. Recently, an extensive 
study has been made by Lagerstrom, Cole, and Trilling (reference l) ; a 
linearized form of the Navier-Stokes equation in which the transverse 
viscous effects and the longitudinal viscous effects can be superposed 
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was integrated. The present paper considers a source or sink flow of a 
viscous, heat-conducting, compressible fluid. This problem is a non- 
linear flow case in more than one dimension for which the viscous effects 
are solely longitudinal. The relative simplicity of source flow is due 
to the fact that although the flow is two- or three-dimensional, the flow 
can be represented by single -parameter equations because of its cylin- 
drical or spherical symmetry and, thus, it has no place for transverse 
viscous effects. The present solution may be considered a generalization 
of one-dimensional, viscous, heat-conducting shock flow to two or three 
dimensions; in such cases, expansion flows as well as compression flows 
are possible. 

In the solution of the Navier-Stokes equation given herein, use is 
made of the second viscosity coefficient. The fundamental developments 
concerned with the . second viscosity coefficient are given in Busemann's 
"Gasdynamik" (reference 2) published in 1931, derived by means of a 
gas-dynamic approach, and in Tisza's paper (reference 3), published 
in 19^2, where the derivation is for ultrasonics and its implications 
for gas dynamics is only briefly mentioned. Since the development given 
by Busemann is very brief and Tisza's development is derived with a 
different approach, a more detailed account of the second viscosity 
coefficient in gas dynamics is given in the present paper. 


SYMBOLS 


X 

distance from source (or 

sink) 


space coordinates normal 

to x 

u 

velocity in direction of 

source (or sink) 

v,w 

velocities normal to u 


V 

— u 

u = 

u 

max 

velocity vector 


T 

viscous stress 


TX 

total stress 



surface of sphere or circumference of circle; also, cross- 
sectional area of tube with slightly varying cross section 


F 
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p 

t 

T 

P 

M 

s 

h 

R 

P 

P ! 

k 

c 

1 

C P 

c v 

7 

Cl 


C = u 
pux 


speed of sound 

density 

time 

absolute temperature 
pressure 

Mach number (u/a) 
entropy 
enthalpy 
gas constant 

first viscosity coefficient 
second viscosity coefficient 
heat conductivity coefficient 
mean molecular velocity 
mean free path 

specific heat at constant pressure 
specific heat at constant volume 
ratio of specific heats ( c p/ c v) 
constant 


0. = a log * 

d log u 


Subscripts: 

o stagnation conditions 

max maximum 


A bar above a quantity denotes a vector except in the case of u 
where u represents u/^. 
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FUNDAMENTAL EQUATIONS AND THEIR INTEGRATION 


The general flow equations are (references 2 and 4): 
The continuity equation is 

div( pv) = 0 


( 1 ) 


The equation of motion is 


Dv . 

p = -dxv n 


( 2 ) 


where it is the total stress tensor which has the following relation to 
the pressure p and the viscous stress tensor t : 


n xx n xy n xz ' 


n y* n yy V 

i5T zx n zy ^zz/ 


/p 0 0\ / T xx T xy T xz 

0 P 0 I + ( T yx Tyy Ty Z 

\0 0 p/ \ T ZX T Z y T ZZ/ 


The energy equation is conveniently written in the form 


div 



(k grad T - v.t) 


= 0 


(3) 


where k grad T represents the heat flow, v * t represents the flow of 
work of the stress tensor t, and the quantity h is the enthalpy of 
the flow. The inner or contracted multiplication v-t (or vfr^) (Tor 
example, reference 5) yields a vector with the components 


(v.t) x 

“ UT xx 

+ 

VT yx 

+ 

WT ZX 

( v* t) , r 

= UT wr 

+ 

VT,„, 

+ 

VT Tr 

'y 

xy 


yy 


zy 

(v* t) 

= UT 

+ 

VT „ 

+ 

V T 

v y z 

xz 


yz 


ZZ 
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The expression div(vr) can he written in the form 




This expression is of the same form as the equation given in reference 4 
with the exception that in the case of reference 4 the stresses include 
the pressure term. 

Due to the symmetry of source flow, instead of transforming equa- 
tions (l), ( 2 ), and (3) into polar form, x and u can he introduced 
for the polar coordinate and polar velocity, respectively. The con- 
tinuity equation hecomes 


In order to express the equations of motion and of energy for source 
flow, the components of the stress tensor are needed (references 2 and 4 ) 
and are 


puF = Constant = C 3 


( 4 ) 





With the use of the symmetry relations for source flow 
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and conditions for irrotationality, the equation of motion becomes for 
three-dimensional source flow 


pu 


du _ d£ du 2 / , _ )(du 2 u\j ^ d /u\ 

" " a£ dx|_^ dx + 3 ^ ^Vdx x/J dxVxy 


and for two-dimensional source flow 


pu 


du _ d£ + aL du 2 ( _ )/du + u\l + AM 

dx ' ' dx dx|f U dx 3^ Vdx x/J dxW 


The equation of energy for three-dimensional source flow is 


c, J-fei + h) - F 

3 dx\2 / dx 


— + 2Pr —£—) + |(p' - P)\ dx 
Cp\dx dx / 3 \ dx 


x 


( 5 ) 


( 6 ) 


= 0 


(T) 


u 2 

For two-dimensional source flow the term u^/x replaces 2 — in the 


energy equation. The Prandtl number Pr is given by pcp/k 

By putting M-' = M- and Pr = ^ in equation ( 
of constant total flow energy is obtained) that is, 


2 

— + h = Constant 
2 

A detailed discussion of the relation between M-' and p is given in 
appendix A. In terms of the energy equation in the general form (equa- 
tion (3)), the constancy of total flow energy is given by 


k grad T - v*t = 0 


The compensation of heat flow and work flow is, of course, in line with 
the statement that no heat is added to the flow. The differences between 
source flow with constant total flow energy and with heat addition are 
discussed in appendix B. 


With the use of the requirement of p' - P for constant total 
energy flow, the equation of motion for three-dimensional source flow is 
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pu 


du 

dx 


- & + 
dx 



du\ 

dx/ 


+ 4p 


d /u 


dx\x, 


and for two-dimensional source flow 


( 8 ) 


pu = _ JLfat + 2P -f(ii) ( 9 ) 

dx dx dx y dx/ dx\x/ 

In terms of the area F, which for two- and three-dimensional flow is 
given by 2nx and 4 jtx 2 , respectively, equations (8) and (9) become 
identical: 


dp 

dx 




do) 


Now, the following substitutions are made in equations ( 8 ) and (9): 



or 


P = P ~ u raax 2 ( 1 - u 2 ) 

where a 2 = 2 _I_i u^d - u^) and is used as the reference 

value. With the assumption that d varies in the form 


d 

Po 



and with the substitution 


T 


7 - 1 
2/R 


u max 2 ( 1 


S 2 ) 


the following equation results: 

= u max C l\f~~ ^ 
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where 


C 1 


I 7 - 1 
i/^V 27R 


Equations (8) and (9) are now transformed to make the coordinate x 


nondimensional by first writing 


djl 

dx 2 


in the form 


d 2 u _ _1_ d^ 1 du 

dx 2 x 2 d log x 2 x dx 


obtained through development of 


d 2 ^ 


d log x c 




+ x“- 



For the transformation from x to log x, multiply equation (8) and (9) 
by x/pj the following equations result: 


For cone flow 


I d lQ g P h . n 2 ) + i u 


1 — du 


2 7 d log x 


+ 200^(1 - u 2 ) 2 (-^ 


TT \2 


7 “ d log x ' “2“ “ \d log x 


lC 2 u 2 (l " u 2 )^ 


2 - 


2C 2 u(l - 52) 2 d -H — 
d log x c 


1 

2 




d log x 


(ID 


and for wedge flow 


7-1 d lop; P (i 
27 d log x 


- n 2 ) + i u — ^ + 2C 2 U 2 (1 - u 2 

7 d log x 


) 2 (_du V 

\d log x/ 


— 2\2 d 2 u 


2C 2 u 2 (l - u 2 ) 2 - 2C 2 u( 1 - u 2 )' 

d log x 


= 0 


( 12 ) 
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The quantity Cg is given by 


c c i JEH. _a_ = ^ Jy - 1 1 

pux [//RTq |f 2 pux a o r 2 pux 

Since u is actually a short notation for u/ u max where u max was 
assumed to be the reference value and since 


^ax _ ./ 2 
a o = 


C2 may be written in the form 


C - M ' 0 ^ax i /y" - 1 _ ^o 
^2 " pux a 0 \j 2 pux 

Now, the mass flow for two- and three-dimensional source or sink flow 
(wedge or cone flow) is given by 


puF = pux2Jt 


and 


puF = pux 2 lw 

For wedge flow the dimensionless parameter Cg> which has the form of an 
inverse Reynolds number, is a constant) whereas for cone flow 


C 




- Constant 


(x) 


or C 2 is proportional to x, the distance from the source or sink. 

For three-dimensional flow no advantage has been obtained by the introduc- 
tion of log x since a dependency on x is hidden in In 

present paper, the simpler case of wedge flow is discussed, since its 
simplicity permits a quick grasp of "the fundamental nature of the flow 
problem. 
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Equation (12) can be simplified by introducing the continuity 
equation in the form 


d log p _ 1 du 
d log x u d log x 


Since in the equation the independent variable log x appears in much 
simpler form than the dependent variable u> equation (12) can be further 
simplified by an interchange of the dependent and independent variables. 

This interchange is achieved by multiplying equation (12) by (^ - 1 °£ — j 

, . d u . 

and expressing by 


d log x 2 


d 2 u 

d / du \ 

_ _d_ 

(d log x ' 

-1 

du 

d log x 2 

d log x\d log xy 

du 

_\ du ) 

1 

d log x 


d 2 log x 

/ du 

du 2 

Id log x 


There fore , 


S'- 


7 1 (1 - u2) + 2C 2 u 2 (l - u 2 ) 2 


+ I 


.2 


7 - 1 
27 


(1 


- u 2 )= + \ ffl 

U 7 J 


1 1 

2 . rr 2x2 _ 


4'2C 2 u 2 (1 - u 2 ) + |"2C 2 u(1 - U 2 ) = 0 


(13) 


or 


p = 


u 


—2 

u 


I* 


liji 2 - 1 

k7C 2 u 2 l/l - H 2 


(7 - 1 V 


V+7C 2 

5 ) 


7-1 


(U0 
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where 


*• = ^- X °g x 

du 

To avoid odd powers in the coefficients for convenience, reduce equa- 
tion (l4) by introduction of a new variable 


j, d0(log u) = d log x 
v ~ d log u _ d log u 


such that it becomes 


0 " = 0 ’ 1 - 


u 


1 - U 2 / 


- 0 


7 + l_2 

- 7 - 1 

,2 y - 1 


4/C; 


7-1 


u 2 /! - u 


772 


0 ' 3 £ 


^7C; 


/l - u 2 
IT 2 


- 1 


( 15 ) 


In the present form, u is always raised to the second power and, thus, 
since the second power term is always positive, the sign of C 2 will 
determine whether the flow will pass in the positive or negative direc- 
tion through a given wedge. Specifically, since 



C 2 will be positive if u and x are of the same sign and negative 
if they are of opposite sign. Since, as previously mentioned, x is 
the distance from the apex of the wedge, a positive C 2 signifies a 
flow through a diverging wedge source, whereas a negative C 2 indicates 
a flow through a converging wedge sink. 

The physical significance of niay be seen by making the substitu- 

tions (see, for example, reference 6) 

fi = 0.b99pcl 


and 


c 


n 7 
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and by applying relations previously used. The relation for C 2 becomes 


c 2 


= K-, 


l fSo 
x u 


where 


k i = o-^i/d 

As may be seen from flow tables (for example, reference 7)> the 
ratio a Q /u will vary comparatively little for a large range of super- 
sonic and even subsonic Mach numbers. For the supersonic Mach number 
range, a Q /u will be represented to a good approximation by a constant 
and C 2 will be proportional to l/x with a proportionality constant 
not far from one. The significance of C 2 , a reciprocal Reynolds number, 
may perhaps be more conveniently seen by thinking of x as the width of 
a given wedge section through which the flow passes. Since the density 
is inversely proportional to the mean free path (reference 6), C 2 is 
about inversely proportional to both the density of the flow and the 
width of flow section. If C 2 = 1 is assumed as a criterion for the 
case where viscous effects due to compressibility are of importance, 
the isolated longitudinal viscous effects will be negligible for hyper- 
sonic tunnels operating in the range of atmospheric stagnation conditions 
even though the minimum section of the tunnel may be small based on 
engineering standards. If the minimum section of the hypersonic tunnel 
is to influence greatly the isolated longitudinal viscous effects with- 
out becoming impracticably small, the mean free path has to be increased 
to the order of the width of the minimum section} that is, the stagnation 
density of the flow has to be correspondingly small. For three-dimensional 
flow C 2 increases with the distance from the source and, thus, the 
viscous effects, will be correspondingly influenced. 


Since the nature of equation (15) rules out a closed integration, 
the integration is performed numerically , which is done conveniently by 
a modified isocline method. The modification consists in requiring 


jjT = Constant along the modified isoclines in the 0',u plane instead 

of requiring 0" = Constant as would be necessary for true isoclines} 
0"/0' represents reciprocal subtangents. The equation of the modified 
isoclines is given by 


0 


7 + 1—p 

7 - 1 U - 1 


(1 - H2) . C 2 

7-1 

. - u 2 


/r 7 + 1 u 2 - 1 1 

/ A 

2 

4(1 - 2TT 2 ) 

± 0" 4 

l jjl - u2) - C 2 u 2 ^ - u 2 

(1 - n 2 )| 

^7 - 1 ^1 - u 2 

1 P 7 7 - 1 Vl - u 2 , 

\ k7 c 2“ 2 ) 

1 w u 2 


(16) 

The integration curves in the 0’^u plane are obtained by connecting the slopes on the isoclines 
corresponding to various values of 0"/0 ! (see figs. 1 and 2). The integration curves in 
0 r ,u plane are most conveniently integrated in the u, x plane step by step by means of the 
expression 


A0 = 0 f A log u 


Since the curves 0(log u) represent an integration of the curves 0 T (log u), they are deter- 
mined except for an arbitrary constant. The arbitrary constant means that the curves 0(log u) 
in figure 3 may be shifted along the x-axis } in other words^ the two -parameter family of curves 
obtained from the second-order differential equation has been converted into a single-parameter 
family and a translation family. This shifting of the solutions along the x-axis could have been 
predicted from the fact that for two-dimensional flow is independent of the distance from 
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the source) for the simple isentropic case, of course, the solutions can 
also he shifted. It can he expected that for three-dimensional source 
flow, a simple introduction of a translation family would no longer he 
possible) thus, a second-order differential equation may have to he 
solved numerically instead of reducing the problem to the much simpler 
solution of a first-order differential equation. 

The integration for each particular (integral) curve was performed 
hy starting at some point of the flow field and following the structure 
of the field from there on. This step-hy-step approach has the aspect 
of an initial-value problem. Since the integral curves in the physical 
plane can be shifted along the x-axis, the curves can he placed in a 
position such that they can he compared with the integral curves for the 
case where the viscous effects due to source flow are zero. In the 
numerical integration a choice had to he made for the magnitude of the 
parameter C2- The values of Cq = -0.1 were selected, since trial 
computations for various values of Cq indicated that in this case the 
longitudinal viscous effects due to compressibility would be compara- 
tively small. Such solutions should he of special interest since they 
can he expected to have certain aspects in common with both the cases 
of zero viscous effects C2 = 0 and finite viscous effects. In order 
to understand fully the contribution of the various aspects, the struc- 
ture of equations (15) or (l6) is analyzed for the neighborhood of C2 = 0. 


DISCUSSION OF SOLUTIONS 

The Neighborhood of C 2 = 0 with a Discussion of C2 = 0 


To investigate the neighborhood of C 2 = 0, equation (15) is multi- 
plied hy C 2 J then C2 appears as a coefficient of 0", the highest- 
order term of the differential equation. The problem becomes one of 
deciding what happens when the coefficient of the highest-order term 
of the differential equation approaches 0. A similar problem arises in 
connection with the development of the bound ary -layer equations and 
equations of shocks in constant cross section where p (instead of C2) 
approaches 0 (references 1, 8, and 9) • (The term shock is used herein 
for compressions with finite p and k as well as for compressions 
with p — >0 and k — »0.) Generally, for this type of problem, the effects 
of viscosity and heat conduction are of such a nature that they smooth 
out the discontinuities existent when the coefficient of the highest-order 
term of the differential equation is set equal to 0. This smoothing of 
the discontinuities applies, of course, only to the immediate neighbor- 
hood of C 2 = 0. For the present case, when C2 = ±0.1, values which 
are finite but small, such a process will be true only approximately, but 
it is very useful in under standing the flow structure. 
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Cg = 0.- The discontinuities for Cg = 0, although obtainable from 

physically plausible considerations, are obtained from equation ( 15 ). 
Equation (15) is rewritten after multiplying it by C 2 



0* 3 


7-1 

47 


\jl~^ 


i[2 





0 


(17) 


If Co is set equal to 0 and 0" remains finite and arbitrary, equa- 
tion tl7) reduces to 


1 

ii2 r - 1 


47 

7-1 


u 2 ^l - u 2 


0.3 7 - 1 \/l - u 2 _ 

47 ° 


u 2 


This equation has one double root 


( 18 ) 


0- 


I, II 


= 0 


and a third root 

? t 1 u 2 - 1 

0 * = 21 ^ 

1 - u 2 

which is the equation for the isentropic curve. The isentrope and the 
0' =0 axis will be the loci of finite and arbitrary values of 0". This 
fact is illustrated in figure 4(a) representing the 0',u plane by fans 
of lines distributed on these curves. The solutions of 0' = 0, or the 

equivalent of — 1°S . . u = 00 ^ actually represent the previously mentioned 
d log x 

discontinuities existent when the coefficient of the highest-order term 
of the differential equation is set equal to 0. For the smoothing of 
the discontinuities by viscosity, the highest-order term of the differ- 
ential equation has to be kept finite which can be accomplished by 
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making 0" equal to °° (see dashed lines in fig. 4(a)). Then, since 
C 2 0" is indeterminate 0' can assume any value to compensate for the 
indeterminate term. Because 0' is arbitrary the connections between 
the isentrope and the discontinuities 0' =0 can he established at 
every point of the isentrope. These connections represent the smoothing 
effects in their embryonic stage. The letters in figure 4 are given to 
correlate the (a) and (b) parts. 

Since the discontinuities do not follow the isentropic law, entropy 
variations occur. (The entropy balance is given in appendix C.) The 
entropy variations for the discontinuities can be understood from the 
fact that the viscous-stress terms and the heat-conduction terms (since 
the total flow energy is constant) are essentially given by 


du 


du 


“ to * C 2 PUX to 


For given values of p, u, and du the viscous effects are constant if 
C 2 and dx/x decrease proportionately down to 0. For an infinite 
value of x (from a source of infinite strength), 0^=0 means flow 
through a constant cross section) the reduction of the source-flow 
equations to shock flow in a constant cross section is given in 
appendix D. 

Since the velocity of source flow depends only on the distance from 
the source, flow discontinuities may also arise because of the termina- 
tion of source flow by an adjoining non-source flow. Depending on whether 
the cross-sectional variation at the junction of source flow and non- 
source flow is discontinuous in the first derivative or is discontinuous 
itself, the terminating flow discontinuities are represented by ^ = 00 

or by velocity jumps. Naturally the source flow does not have to be 
terminated by discontinuities in the cross-sectional variation at the 
junction or its first derivative) rather, termination of source flows is 
possible for discontinuities in any derivative. At the junction of non- 
source-flow boundaries with source-flow boundaries, the symmetry of 
source flow no longer holds. Since the present calculations deal with 
pure source flow, they neglect the disturbance of the flow symmetry due 
to the termination of source flow. Such unsymmetrical effects will tend 
to disappear for slightly diverging source flows (or slightly converging 
sink flows) since for them the flow velocity tends to be constant in a 
given flow cross section. 
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The discontinuity = °° at sonic velocity is due to the fact 

that a nonviscous source or sink flow is also terminated or limited as 
it cannot pass beyond the location where the sonic velocity occurs. For 
two-dimensional source flow the so-called "limiting line” (reference 10) 
of the flow is a circle at sonic velocity. Because of the use of sym- 
metrical conditions of source flow, in the present case the limiting 
line appears as a point of the isentrope at M = 1 in the physical 
plane (u, x) . 


The immediate neighborhood of C2 = 0.- In investigating the 

immediate neighborhood of C2 = 0 (see figs. 5 and 6), the purpose is 
not to give a detailed mathematical account of the Neighborhood problem” 
but rather to inspect equation (15) for certain features which should 
give an insight into its solutions without requiring actual integration. 
Since C 2 is now finite, though in the immediate neighborhood of 0, the 
right side of equation (17) will no longer be 0 for finite and arbitrary 
values of 0” and thus some of the simplicity inherent in the solution 
for C 2 = 0 will be lost. Equation (15) indicates, however, that for 
finite values of C 2 the equality of 0 of its right side (which also 
results in 0” = 0) will give certain useful relations. Namely, the 

0” = 0 axis will not only be a locus of constant slopes — $$ * _ = 0” 

d log u 

but also an integration curve (in the 0 f ,u plane) since its own slope 

is a YQg " = 0* Since the integration curves of equation (15) cannot 

cross each other, as the coefficients of equation (15) are single-valued 
functions of u, the fact that the integration curve 0* = 0 is a 
straight line can be used for orientation among the infinite group of 
integration curves in the 0 , ,u plane. Through inspection of equa- 
tion (15), orientation is also possible in the direction field given by 


the curves along which _ = 0” = Constant (in the 0',u plane). 

Since the direction field given by equation (15) is continuous, each 
curve of the system of all curves of 0” = Constant furnishes a dividing 
line for the slopes 0" in the plane 0 ? ,u. Similar to what was done 
previously for the integration curves, one particular dividing line for 
the slopes is chosen for orientation among the lines of constant slope 0” 
in the 0 T ,u plane. The dividing line 0” =0 is chosen since in that 
case equation (15) yields as one solution the 0* = 0 axis. The general 
trend of the other two solutions of equation (15) for 0” = 0 Cor of the 


slightly modified equation (l6) for 
hood of C 2 = 0 
are given for C 


£-°) 


for the immediate neighbor- 
= 0 curve s 


is taken from figures 1 and 2 where the 
= ±0.1. Further, for orientation in the direction 
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field, the fact contained in equation (15) can be used that as 0' 
tends to infinity, 0" (the slope in the 0',u plane) will do the same. 
The only exception is that when the coefficient of the term 0'J, which 
has the highest degree in equation ( 15 )> becomes 0, then 0" will be 
indefinite for an infinite 0' . 

The direction fields in the 0',u plane for the positive neighbor- 
hood C 2 = +0 (fig. 5) and the negative neighborhood C 2 = -0 (fig. 6) 
of C 2 = 0 may now be constructed. In order to obtain a clear illustra- 
tion, the size of the regions of finite slopes (furnishing the smoothing 
effects) has been exaggerated in figures 5 and 6. The following features 
in the neighborhood solutions may be singled out. For the immediate 
neighborhood of C 2 = 0 integration curves exist which essentially 

remain in the neighborhood of the isentropic curve for C 2 = 0. The 
direction fields also indicate that integration curves due to compres- 
sion shocks or due to flow changes caused by termination of source or 
sink flow will practically join the neighborhood solutions of the 
isentropic curves. In figures 7 and 8 a summary of the behavior of the 
integral curves for the immediate neighborhood of C 2 = 0 is made. The 
arrows in these figures indicate the velocity increase and decrease corre- 
sponding to source (C 2 = +0) and sink (C 2 = -0) flow, respectively. 

The curves of figures 5 to 8 are numbered (as far as possible) in 
accordance with the comparative curves in figures 1, 2, and 3. Curves 1, 
2, and 3 apply to the case of C 2 = +0 or source flow. Curve 1 repre- 
sents an expansion for infinite source flow and is essentially due to 
the smoothing by viscosity between the discontinuity limiting the flow 
at sonic velocity and the supersonic branch of the isentrope. Curve 2 
represents a compression shock in infinite source flow. Behind the 
shock, subsonic velocity exists which may be decelerated to zero in the 
diverging source flow. Curve 3 represents the smoothing by viscosity of 
expansion flow due to termination of source flow. The curves 4, 5> an< i 6 
represent corresponding effects for Cp = -0 or sink flow. The curves 
are discussed in greater detail for flows with C 2 = fO.l, presented in 

figures 1, 2, and 3* 

Now, a comparison can be made between the flow structures for the 
immediate neighborhood of C 2 = 0 and for C 2 = 0. For the immediate 
neighborhood C 2 is finite, though small, and thus on a strict math- 
ematical basis the direction fields of the differential equation are such 
that all the solutions are represented by separate curves. Figures 5 
and 6 for the direction fields of the immediate neighborhood, however, 
indicate that solutions due to shocks or termination of source-or-sink 
flow boundaries will join for all practical purposes the neighborhood 
solutions of the isentropes. For C 2 = 0, the junctions between discon- 
tinuous changes and the isentrope are correct on a strict mathematical 
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basis. For all practical purposes the neighborhood case still has 
certain aspects which were exact for the case of C 2 = 0. The solutions 

for C 2 = ±0.1 given in figures 1, 2, and 3 indicate that the viscous 
smoothing effect of discontinuities, due to termination of source or 
sink flow o;r due to shocks (p — k — > 0 ) in infinite source or sink flow 
(curves 2, 3, 5, 6, and 7), will approach the expansion curves 1 and 4 
for infinite source or sink flow with sufficient rapidity so that they 
may be joined with them for all practical purposes. In other words, the 
plotted solutions for the finite values C 2 = ±0.1 are close enough to 
those for C 2 = 0 to show the typical aspects of the immediate neighbor- 
hood solutions C 2 = +0 and C 2 = -0. 


The numerical integration, as previously stated, was performed by 
starting at some point of the flow field and by following the structure 
of the field from there on. In integrating this initial-value problem, 
only two curves in figures 1, 2, and 3, curves 1 and 4 (which may be 
shifted), exist which extend from the neighborhood of the source or sink 
to infinity without being interrupted by curves rapidly deviating essen- 
tially in an exponential manner. For finite values of C 2 the rapidly 
deviating curves have to be actually separate curves from 1 and 4, but 
for the small values of = ±0.1 a junction will still be correct for 
all practical purposes. Since the exponentially deviating curves approxi- 
mately join the more gradual curves 1 and 4, integration is advisable in 
the direction opposite to that in which the exponential deviations are 
expected to occur. Otherwise, difficulties will be encountered in fitting 
these solutions into one particular exponential deviation (with certain 
initial values 0 and 0 ? ) of the gradual solutions as exhibited by 
curves 1 and 4. Similar difficulties were found in the analysis of refer- 
ence 11, which deals with one -dimensional shock flow for arbitrary Prandtl 
numbers) other fundamental aspects of that flow structure are presented 
in a very clear manner in reference 12. 

Curve 1 in figures 1 and 3 is the equivalent of curve 1 in figure 5* 

It thus essentially represents the viscous counterpart of isentropic expan- 
sion flow for the case of an infinite diverging wedge. Curve 1 was shifted 
along the x-axis into such a position as to make possible a comparison 
with isentropic flow through the same area ratio. fit should be noted 


The Case of C 2 = ±0.1 


that log x^ - log = log 
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Curve 2 represents a compression shock for infinite source flowj 
its velocity distribution is essentially the same as in the well-known 
shock in constant cross section. Curve 3 represents an expansion due to 
termination of source flow. Curve 4 represents the expansion flow 
through an infinite converging wedge. The flow starts at zero velocity 
and goes through the subsonic range and then essentially smooths out the 
region which is abruptly limited, at sonic velocity for isentropic flow. 
The surprising thing about this expansion flow is that it penetrates 
beyond the sonic section or the minimum section for isentropic flow. 

The meaning of this penetration is discussed subsequently. Curve 5 
represents a compression shock in converging supersonic infinite sink 
flow. Such a shock is actually unstable, but such considerations do 
not enter into the present problem. Curve 7 represents another compres- 
sion shock for infinite sink flow which is unusual in that it is entirely 
in the supersonic region. Its existence is possible because expansion 
curve 4, for converging-wedge flow, penetrates into the supersonic region. 
Finally, curve 6 represents an expansion due to termination of sink flow 
in the subsonic-flow region. Shock 7 follows more closely the trend of 
supersonic isentropic compression than does shock 5* Curves which will 
follow this trend even more closely could have been given. In connection 
with the possibility of existence of such various compression curves, it 
should be noted that for nonviscous flow through converging wedges, 
solutions do exist with shocks and without shocks. 

Since the structure of the flow field has been treated, certain 
general aspects can now be discussed. For example, in figure 3 the 

curves representing viscous flow solutions have a smaller slope than the 

corresponding curves for C 2 = 0. This smaller slope is immediately 
apparent for curves which for C2 = 0 are represented by disconti- 
nuities — = <» or lines of discontinuities. Such is the case for the 

dx 

shocks, curves 2 and 5, and for source- flow terminations with the appro- 
priate boundary changes. (See curves 3 and 6.) In order to indicate 
the fact that curves 1 and 4 have a smaller slope than their isentropic 
(C 2 = 0) counterpart, the isentrope has to be shifted along the x-axis. 

The smaller slope of the viscous-flow curves may also be expressed by 
the statement that the viscous flow requires a larger area ratio to go 
through a given velocity change than does the corresponding flow for 
C 2 = 0. This behavior could have been more or less expected except for 
the case of viscous sink flow converging to the section where M = 1, 
given by curve 4. Namely, viscous effects alone are known (for example, 
reference 2) to cause a given mass flow approaching M = 1 to require 
a larger cross section than the same isentropic mass flow. In other 
words, for viscous effects alone the flow approaching M = 1 can only 
pass through smaller area ratios than the same isentropic mass flow. 

The reason for the fact that the flow given by curve 4 requires a larger 
area ratio is based on the fact that heat can be conducted downstream in 
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the present case. In summary, this behavior means that the familiar 
heat- conduction effects combined with the peculiar longitudinal viscous 
effects alone (no transverse viscous effects) may cause a given mass 
flow to pass through a cross section smaller than the isentropic minimum 
section; in other words, a larger mass flow may pass through a given 
isentropic minimum cross section than for isentropic flow. This problem 
is discussed in greater detail in a subsequent section and in appendix C. 

Now, a few general aspects of the connection of the shocks and 
terminating solutions with the solutions for infinite source or sink 
Hows, curves 1 and 4, can be discussed. Figures 1 to 3 show that for 
flow through diverging wedges, source flow, the exponential deviations 
from the expansion curve 1 become larger in the direction of increasing 
velocities (see curves 2 and 3). For converging- wedge flow, sink flow, 
the exponential deviations in the neighborhood of the expansion curve 4 
increase with reduced velocities (see curves 5 and 6). Since both expan- 
sion flows are, of course, in the directions of increasing velocities, 
the eHects of shocks and source or sink terminations will penetrate both 
upstream and downstream. The nature of the criterions for upstream and 
downstream penetration can be seen from the simple example of the shock 
in constant cross section; namely, since only shock from supersonic to 
subsonic velocities are possible, the effects of a compression shock will 
penetrate upstream on the supersonic side and downstream on the subsonic 
side. The same behavior is exhibited for the present case of source 
flow. For example, curves 2 and 3 indicate that the effects of shock or 
termination of source flow penetrate upstream in the supersonic region 
of expansion curve 1 of the source flow; whereas curves 5 and 6 indicate 
that the same effects penetrate downstream in the subsonic region of 
expansion curve 4 of the sink flow. The shocks show, of course, both 
types of penetration. Effects of shock 2 penetrate downstream and join 
the subsonic part of the isentrope. For the present cases of source or 
sink I low, exceptions exist to these criterions of subsonic and super- 
sonic penetration. For example, shock 7 penetrates downstream into the 
supersonic velocity region of curve 4. This particular supersonic 
velocity region is characterized by the fact that curve 4 is continued 
to u values which would have been prohibited for the case that C2 = 0. 
The behavior of the exponential penetration may thus be summarized. In 
the flow regions which were not prohibited in the case of C 2 = 0, flow 
changes (for example, like those of shocks in infinite wedge flow' or due 
to termination of source or sink flow) will penetrate upstream in super- 
sonic regions and downstream in subsonic regions. In regions pro- 
hibited for C 2 = 0, the nature of the penetration effect if given by 

the character of the flow of the corresponding expansion curve before it 
enters the region prohibited for the case of C 2 = 0. 

Certain limitations of the solutions in figures 1, 2, and 3 exist 
which are of importance for the interpretation of these results. For 
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example, curve 1, the expansion curve for infinite source flow, will 
not reach the maximum flow velocity hut rather will he asymptotic to 
a velocity equal to about 0.7 of the maximum velocity. The existence 
of an asymptote for viscous expansion flow through the diverging cross 
section of source flow or the existence of the asymptote itself as a. 
monotone solution could have been expected. However, the asymptote is 
reached for the case that the normal stress ityy = 0, which is a condi- 
tion for which the Navier-Stokes equation reaches its limit of 
applicability. 


The Limits of Applicability of the Navier-Stokes Equation 

Based on the kinetic theory of gases (see discussion of Burnett's 
work in reference 13 (especially p. 271) ) , the Maxwell distribution is 
regarded as a first approximation of the Boltzmann equation for the 
general distribution function and the viscous-stress terms of the Navier- 
Stokes equation are regarded as a second approximation. The limit of 
applicability of the Navier-Stokes equation is considered to be reached 
or passed when the terms of the third approximation reach the same order 
of magnitude as the terms of the second approximation in the Navier- 
Stokes equation. Actually, of course, no sharp limit exists for the 
applicability of the Navier-Stokes equation, but the reaching of the 
same order of second- and third-approximation terms is a sufficiently 
broad criterion or barrier to be generally accepted. A discussion of 
limit of applicability of the Navier-Stokes equation is furthermore 
simplified by the fact that the ratio of the third-approximation terms 
and the second-approximation terms is of the same order as the ratio of 
the terms of the second approximation and the first approximation (see 
also reference 14, p. 453). For the present case, the limit of applica- 
bility of the Navier-Stokes equation is reached when the viscous stresses 
reach the order of magnitude of the fluid pressure. 

Since the present solution is based on the condition P' = b, the 
normal stresses Tt xx and Ti y-y are ( rePeI>ences 2 and 4) : 


n xx p " dx 

and 

"yy = p _ ^x 
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Actually,, a choice exists between two limits of applicability of the 
Navier-Stokes equation; one for 



and the other for 



In terms of the variables used in the present paper^ these two break- 
down criterions are expressed by 



Breakdown 


r - 1 ~ 


and 


I2n £uj 


^/Breakdown 


^ 7( "2 d log x _ ^ 

7 " 1 ]fl - u 2 d log u 


Since /l - u 2 is positive because it is related through the 
temperature T to the necessarily positive viscosity, the first break- 
down criterion p = 2^ has a meaning only for positive values of C 2 , 

that is, for flow through a diverging wedge. For C 2 = 0.1 and 7 = 1.4 
(for air), u =_0.7095- This value is the asymptotic value for which 

the gradient — is 0 for the expansion flow through an infinite 

diverging wedge (curve 1, figs. 1 to 3)- Also, all expansion curves for 
diverging wedge flow will have to be cut off at u = 0.7095. The condi- 
tion p = results also in flyy = 0,* thus, the stress normal to 

the streamlines is 0. The flow of the gas is thus no longer forced to 
follow the diverging streamlines, or, in other words, if the divergence 
of the wedge were to be increased, the flow would no longer follow. 


Two possibilities, those of expansion and compression, exist for 
the second breakdown criterion p = 2 a ||. In table I, the ratios of 


I2p || 


to 


2H 


du 


are given for the various integration 


Plot 


dx /Breakdown 
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curves in figure 3* These ratios are equal to ratios of 


d log 

d log x /Breakdown 


f d log u | 

\d log x /piot 


or 



Plot 


where 0’ is the ordinate in figures 1 and 2. 


Table I indicates, for example, that for compression shocks in fig- 
ures 1, 2, and 3, breakdown of the Navier-Stokes equation will occur in 
certain flow regions within the shocks similarly to the well-known case of 

one -dimensional shocks. fit should be pointed out that the breakdown 
P ' 

criterion =1 is less severe for compressions than for expansions 

a* ^ \ 

dx \ 

for which it coincides with it xx = O.j The table, furthermore, shows 


that for the expansion curve 3 based on source flow, breakdown due to 

du 

p = 2 m-^ (u = 0.7095) occurs before breakdown due to p = 2 m- The 


table also indicates that the breakdown of the infinite source (curve l) 
and sink (curve 4) solutions will occur in those flow regions in which, 
for isentropic flow, no solutions would have existed at all. Such 
behavior was already suggested in previous considerations. 


As is indicated by C 2 = Constant (x) for three-dimensional source 
flow with constant total energy, the limit of applicability of the Navier- 
Stokes equation will be reached sooner than for two-dimensional flow 
(C 2 = Constant) if x increases and later if x decreases. Since the 
equations for two- and three-dimensional flow differ essentially only by 
the dependence of C2 on x, for flow changes in small regions (small 
changes in x) , like those occurring in the regions of shocks and flow 
termination, only a small difference will exist between two- and three- 
dimensional flow. 


Flow through a Curved Minimum Section Joined to a Sink Flow 

Since a compressible source or sink flow cannot reach the source 
or sink, the problem arises concerning its continuation. The problem 
is of special interest in the case of the sink flow given by curve 4 
which is able to penetrate beyond the M = 1 section. In this section 
flow in a flow filament through a curved minimum section joined to the 
sink flow is treated. Since the expansion flow of curve 4 passes to 
smaller cross sections than the isentropic minimum section, it would 
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appear plausible on first thought that the flow would also pass through 
the curved minimum section. In a curved minimum section joined to the 
sink flow or converging wedge flow, the cross section is further reduced, 
however, below its value at the junction and, thus, a separate investi- 
gation has to be made by solving the differential equation for arbitrary 
boundary shape s . 

Since the solution would result in a laborious numerical task, an 
attempt will be made to gain information by a study of the nature of 
the differential equation without actually solving it. For this purpose, 
the two -dimensional -flow equations will not be investigated, but equa- 
tions simplified by assuming the flow to be quasi -one -dimensional with 
slightly variable cross section (filament flow) will be investigated 
instead. (The actual two-dimensional flow through the curved minimum 
passage, unlike the wedge flow, also contains transverse viscous effects. 
These effects are, however, eliminated for the simplified case of 
filament flow.) For this case, the flow equations may be written in 
the same form as for wedge flow, with the only difference that F is 
no longer equal to 2nx but that the variation of F as a function 
of x and the flow variables are arbitrary. The fundamental equations 
are transformed in such a manner that the effect of the flow variables 
on the variation in cross section may be easily recognized. 

After the expression 


is substituted in equation (10), the equation of motion is obtained in 
the following form: 


*XX « p - * ^ 



The continuity equation is 


puF = Constant = 


Also, for the present case of constant -energy flow, as previously 
determined, 


T = W 2 d " * 2 ) 


^ u maxCi/l " u2 
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Now, the equation of motion may be rewritten as 

— t? 1 Z2 d u d log F 1 dlt xx . 


- pu s + ^ u d - i d < y ~ + ■ u2 H5E 


dx UwiQ „2 dx 

u max 


= 0 


or 


d 2 log F 


dF\ 1 du + 1 


1 d,T xx 


1 2C 1 - u 2 djyuF dx - u 2 u u. 


=2 - „ 2 dx 

■max 


= 0 


where 


C ^ = PuF 

The relation between and p = pRT may be written as 


z__i (1 . tf) . -5 sl + § 


^ax 


du _ 1 

dx " i/m? 


47c 


1 ^3(1 “ n xx 1 

uF 


>W 22C 1 


If the preceding expression for — is substituted into the rela 

QX 


tion for - log L^ -, the modified form of the Navier-Stokes equation is 


dx 


obtained: 

d 2 log F 


dx c 


,2Cim 


ur 


2C- 


dF\ 

1 ( 

1 — 1 

1 


/l - u 2 ' 


1 

d v _ 

0 

u 2 
u max 

dx 


^ ic_ 1 


XX 


U, 


•max 2 2 ^1 


( 19 ) 
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For determining the connection of a curved minimum passage to the 
sink flov (given by curve k) , it should be recalled that the reason 
that the sink flow could he continued only up to a certain distance 
beyond the isentropic minimum section was the breakdown of the Navier- 
Stokes equation. The shape of a nonisentropic minimum passage is thus 
limited by the requirement that > 0. Since, beyond the isentropic 

minimum section, the Navier-Stokes equation is close to the breakdown 
criterion (see table I, curve k) , the condition for the minimum passage 

/ x dn xx 

will be introduced in equation (19) that tt xx = Constant or = 0. 


In order to determine in principle if the flow will pass through a 
smaller than isentropic curved minimum section, it appeared more expedient 
to investigate the fundamental behavior of equation (19) by simple inspec- 
tion rather than by integration. Since the entropy balance (appendix C) 

diX 

is indirectly contained in the modified Navier-Stokes equation, — has 

to be positive) furthermore, for the converging part of the minimum 
dF — 

passage ^ is negative. Thus, if u does not exceed its physical 
limits and thus invalidates the result (this problem is discussed subse- 
\ d^log F 

quently), § — will remain positive. In order to understand the 

dx^ 


meaning of this condition, the variation of the cross section 


sponding to the simple cases 


d^log F 


= Constant and 


d^log F 


dx^ 


be briefly discussed. Integration of 


a 2 iog f 


dx^ 


= Constant = A 


dx £ 


F corre- 
= 0 will 

results in 


F 


AX^ -n 

^r+Bx 
= Ce 2 


where A, B, and C are constants. 

Now, since at the junction of the converging wedge and the con- 
necting piece d ^ is negative, B will be negative. The quantity A 

2 

is positive since A F has to be positive. The quantity e raised 

dx 2 

to the negative value of a power of x will be represented by a curve 
asymptotic to the positive x-axis. However, the quantity e raised to 
the positive value of the power of x will be asymptotic to the negative 
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x-axis and, thus, will positively increase. Since the positive contribu- 
tion is connected with x 2 and the negative contribution with x, the 

positive contribution will be the deciding one and d . P-P i a -H will form a 


dx c 

variation of cross section which has a minimum section and increases 

downstream of it. The condition d - ■ 1 - ° - | ^ = 0 which does not agree with 

du . 

the requirement of a positive — will yield 

dx 


F = Ce Bx 


where a negative B represents a curve asymptotic to the positive 
x-axis . 

Finally, it is investigated by inspection whether u exceeds its 
physical limits. A lack of a limiting value of u, coupled with the 
continuity condition and other available relations, will indicate that 
the other flow variables will not reach limiting values. Thus, the 

problem arises whether the requirement of a positive may cause u 

to increase to its limiting value u = 1. If u approaches 1, the 

first term in the expression for upon eliminating the brackets, will 

go to 0 and the second term will go to infinity} since the second term is 
subtracted from the first term, this behavior means that as u goes to 1, 

n— 

— becomes negative, which again indicates that actually u = 1 will not 

be reached. A check of whether ^ will actually remain positive as 

dx 

indicated by the entropy balance (indirectly contained in equation (19)) 
can be made to some extent by simple inspection of equation (19)- (Such 
a check will also show if any unexpected irregularities occur in equa- 
tion (19).) The term in parentheses 


( 7 - 1 — 1 - u^ Const ant ^ 

57 ~ 3 T3F " 2Cp J 

indicates that for u = Constant and decreasing F the first term in 
the parentheses will increase, and since the second term is a constant, 

du will remain positive. The fact that the first term includes 1 - u 2 
dx _ 

does not affect the considerations since u = Constant. These considera- 
tions do not yield information concerning the minimum section itself for 
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dJT 

which — =0 or for increasing cross sections. For these conditions, 
dx 


as for all other states at cross sections smaller than the isentropic 
minimum section, the entropy balance, as previously stated, yields 
sufficient information for the present purpose. For increasing cross 
section, the first parenthetic term on the right side of equation ( 19 ) 
will tend to reach 0 after a while. A more exact investigation of the 
case of variable cross section is not within the scope of the present 
paper. 

Similar to the flow passage through a minimum section of sink flow 
alone, for the case of flow through a curved minimum section joined to 
a sink flow, the flow passage through a cross section smaller than the 
isentropic minimum cross section is possible. In other words, a larger 
mass flow may pass through a given isentropic minimum cross section than 
for isentropic flow. The significance of the result can be amplified 
by a comparison of the conditions of mass flow through a minimum 
section for the case of free -molecular flow and the present viscous - 
flow case for which the transverse viscous effects are neglected. Since 
the results obtained in the paper apply not only to minimum sections 
of source or sink flow but to curved minimum sections in general, an 
especially simple case of free -molecular flow can be chosen for a 
comparison. It is the case of effusive molecular flow from a vessel 
(reference 6 ), for which, similar to the present viscous flow, no shear 
at the flow boundaries will exist. The reason that free -molecular 
source or sink flow could not be chosen for comparison is that the 
conditions of molecular chaos on which the Maxwell distribution is based 
are no longer fulfilled for source or sink flow where all molecules have 
to leave or enter a given point (three dimensional) or line (two dimen- 
sional) . The mass flow for the special case of effusive free -molecular 
flow is (reference 6 ) 


where the subscript o refers to stagnation conditions. The isentropic 
mass flow through a minimum section is given by (see reference 7 ) 


Thus, the "very viscous" effusive free -molecular flow has a smaller mass 
flow than the isentropic flow. In contrast, the present solution of the 
Navier-Stokes equation, dealing only with the longitudinal viscous effects, 
may have a higher mass flow. Thus in this respect, the viscous flow 



(for air, 7 = lA) 


°.579a 0 p 0 


(for air, 7 = 1.4) 
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does not follow the trend of the free -molecular flow. The main reason 
for this exceptional behavior of the flow apparently lies in its ability 
to send energy ahead (by heat conduction) , which cannot occur for 
isentropic flow. The quantitative nature of this trend still has to be 
determined, namely, whether the trend towards larger mass flow indicated 
by the present solution of the Navier-Stokes equation will yield much 
higher, or negligibly higher, maximum values than the isentropic flow 
before the mass flow drops to the value for free -molecular flow. In 
order to check this trend quantitatively, a more exact relation than 
the Navier-Stokes equation should be investigated. 


CONCLUSIONS 


A study of the solution of the Navier-Stokes equations for source 
and sink flows of a viscous, heat-conducting, compressible fluid yields 
the following conclusions: 

1. The fundamental effect of viscosity and heat conduction is a 
smoothing out of the flow discontinuities which exist when the viscous 
effects due to source flow are zero. Such smoothing effects of flow 
discontinuities are well-known in boundary -layer theory and one- 
dimensional shock flow theory. 

2. The influence of the familiar heat-conduction effects combined 
with the longitudinal viscous effects alone (no transverse viscous effects) 
on a flow with constant total energy may cause a larger mass flow to pass 
through a given isentropic minimum cross section of a sink flow than for 
isentropic expansion flow. The reason for such a surprising behavior is 
that heat may be conducted downstream. The same effect applies also to 
flow through a curved minimum section joined to a sink flow. This trend 
does not follow the comparable case of free -molecular flow through a 
minimum section, since the free -molecular flow will have a smaller mass 
flow than the corresponding isentropic case. The quantitative nature of 
this different trend should be investigated with relations more exact 
than the Navier-Stokes equation. 

3 . The solutions for viscous, heat-conducting, compressible fluid 
require larger area ratios for the passage of a fluid through a given 
velocity change than do the corresponding solutions for the cases that 
the viscous effects due to source or sink flow are zero. 

4 . The order of magnitude of the longitudinal viscous effects is 
indicated by a dimensionless parameter having the form of a reciprocal 
Reynolds number. This parameter is independent of position from the 
source or sink for two-dimensional flow and increases with the distance 
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for three-dimensional flow. The parameter indicates that the isolated 
longitudinal viscous effects should be negligible for hypersonic tunnels 
operating in the range of atmospheric stagnation conditions even though 
the minimum section of the tunnel may be small based on engineering 
standards. If the minimum section of the tunnel is to influence greatly 
the isolated longitudinal viscous effects without becoming impracticably 
small, the mean free path has to be increased to reach the order of the 
width of the minimum section; that is, the stagnation density of the 
flow has to be correspondingly low. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., October 2 6 , 1951 
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APPENDIX A 

DISCUSSION OF THE SECOND VISCOSITY COEFFICIENT p' 
FROM THE GAS-DYNAMIC APPROACH 


The expression for the second viscosity coefficient p', given 
by Busemann in reference 2, is 



(Al) 


where n represents the number of degrees of freedom and is obtained 
by comparing the expressions for the normal viscous stress t xx based 
on continuum theory and kinetic-gas theory. The two expressions are 


and 

T xx = -2K 2 Pcz[|^ - + ^ + if )] (A3) 


where K 2 is a constant. 

These equations indicate that if only the three translational 
degrees of freedom of a molecule have to be considered, that is, n = 3> 
p' will be 0. Since the existence of p' is due to excitation of 
the internal degrees of freedom, a brief discussion of its relation to 
the well-established time-lag effect should be of interest. The use of 
this effect in gas dynamics is discussed in references 15 and 1 6 , and 
the application to ultrasonics is given in reference 17 . (Also see 
bibliography of reference 17*) The expression of p' through the time- 
lag effect is obtained by modifying the effect in such a manner that it 
is in agreement with the fundamental assumptions inherent in the Navier- 
Stokes equation. Tisza shows (reference 3) that for the case of ultra- 
sonics the time lag may only be incorporated in the Navier-Stokes 
equation if cut « 1, where t is the time lag of the internal degrees 
of freedom and cn is the sound frequency. For gas dynamics, this rela- 
tion is conveniently written in the form 


tflow 
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that is, as the ratio of the time lag of the internal degrees of freedom 
and the time in which the change in flow variables, for example, the 
velocity, occurs. In other words, this relation indicates that tp n t 
may be large, provided "tfiow is much larger. The limit of applica- 
bility of p ' (analagous to the previously discussed limit of applica- 
bility of the Navier-Stokes equation) is assumed to be reached when tj_ n ^- 
is of the same order as ffp ow * 

The detailed derivation of equation (A3) based on kinetic-gas 
theory is presented as follows: Equation (A3) represents the normal 

stress in the x-direction acting on a three-dimensional compressible 
fluid element due to a small deviation from the equilibrium state or 
the Maxwell distribution, for the case that the deviation causes no 
transverse viscous effects. It is convenient to show first the effects 
of a deviation for a gas with one degree of freedom. The equilibrium 
state is given by the pressure (reference 18 which is helpful for an 
understanding of the entire subsequent derivation), 

Px = Px c x 2 

The subscript x indicates the x-direction. Note that p x is a line 
density. The deviation from the equilibrium state is 


dp x = 2p x c x dc x + 


c x^ p x 


For the definition of c x ^ in view of the general nature of the problem 
no distinction is made between various methods of averaging. The expres- 
sion for the deviation dp x is simply made to apply to a three- 
dimensional gas by using the volume density p instead of just the line 
density p x . The average velocity c x is not affected by this general- 
ization to three dimensions, since the three-dimensional Maxwell dis- 
tribution may be obtained by superposing three independent one -dimensional 
equilibrium distributions in the x-, y-, and z -directions . This inde- 
pendence is also expressed by 


c 


2 

x 


c 


2 
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The deviation dp x , which makes use of the correct volume density, 
may be written 


The term dp may be expressed by the continuity equation in the form 


In order to make the results obtained by kinetic theory comparable to 


Since the stress represents only the irreversible effect of 

the deviation from the equilibrium state, the effect of the reversible 
or isentropic deviation has to be subtracted from the preceding expres- 
sion. In contrast for stress based on purely transverse viscous effects 
a deviation from equilibrium can only cause irreversible effects. The 
reversible effect of the deviation is given by 



= -p div v 


dc x 

those of continuum theory > is expressed in terms of the velocity 

's. c x 

gradient For the case of small deviations, the following relation 

is true: 



du 


u 


or since u 


dx 

dt 


c x dx 


Substituting these relations in (hp x )vol 6 dves 



/ bp\ _ dp 
\ P /rev p 
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With the use of the continuity condition and 


7 


2 + n 
n 


the following equation can be written 


( dp) = - ^ + n p div v dt = - -- — pc^div v dt 

rev n n 3 

Subtracting (dp) rev from (^Px)vol as P rev i° us ly suggested gives 
( d Px)vol - < d P>rev " - k pc 2 f^ g; + ( x - ^-^) dlv ”] dt - 

= - 2 (l pc2 )(s - i div 7 ) dt 

In order to obtain the stress t xx , as given by equation (A3), 
integrate the preceding equation with respect to the time t. For the 
present purpose, the integration can be performed to a good approximation 
by changing dt to At which can be done for the following reasons: 

The integration with respect to time indicates that the magnitude of the 
deviation t xx from the reversible reference level will depend on the 
time the deviation can be regarded as independent of the effects of 
collisions. The simplification of the integration by the use of At 
means that a certain average value is assumed for the deviation which 
is abruptly terminated after a time At when the energy due to the 

velocity change in one direction ^see term is suddenly distributed 

over all degrees of freedom ^see term ^ div v j. The time At, in 

accordance with basic concepts of kinetic theory, is approximately 
given by the ratio Z/c where Z is the mean free path. The intro- 
duction of a finite time does not violate the assumption of small 
deviations in the Navier-Stokes equation. As previously shown, the 
immediate neighborhood is a relative concept, based on the requirement 
1 2.nt 

that the ratio « 1. For this case the time lag of the internal 

^f low 
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degrees of freedom equals that of the translational degrees. The 

stress t xx may he written thus 

T xx = - 2 (l - S dlv7 ) < M > 


The expression is identical with equation (A3) if K 2 is set equal 

to 1 / 3 . 

In accordance with a suggestion given in reference 2, the preceding 
equation (a 4) is generalized by taking into account the fact that not 
all the molecular degrees of freedom will be excited after a single 
collision or after a single mean free path: 



where the subscripts 2 and 3 refer to other than translational 
degrees of freedom and P is given by K 2 pc2. Now, an expression for 

I i* can be obtained by comparing the expression (A5) with equation (A2), 
rewritten in the following convenient form: 


T xx = 


35 + Kv" lldlv v 


Comparison of the preceding two expression for t xx yields 


= . 1 + h(i_ . + h(±. . 

3 ^1 \ n l n 2/ Z l\ n 2 n 3, 


lftLl 
3 \h 


1 


* 4 - • . 
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For |-i ' 


til 

P 


3 


p2. 

(1 

l\ 

1 3 1 

(l 

l\ 

\h l 

I 11 ! 

nj 


\ n 2 ' 

— + . . . 
n 3 / J 


or 



For purposes of demonstration it is of interest to make a rough 
estimate of the numerical value of the expression for p*/p for a 
diatomic gas at room temperature. The three translational degrees 
(n-^ = 3) and the two rotational degrees (total iig = 5 ) are excited 


after one collision 



For room temperature, the number of 


collisions or the number of mean free paths necessary to excite the 

vibrational energy is very large, for example, of the order of 10^ (for 
a gas without impurities). The vibrational energy at room temperature 
is, however, only a small fraction of a percent of the vibrational 
energy at equipartition; thus, in terms of the expression for p»/p not 
two vibrational degrees (total = 7) are excited but only a very small 

fraction (total n^ = Very small fraction over 5^* The contribution of 

the vibrational term, the product of a large and a small quantity, can 
thus be expected to be of the order of magnitude of one. Note further 

that because of the limit of « 1 to the validity of p* for 

^flow 

very rapid flow changes like shocks (t flow is small), it may be impos- 
sible to incorporate the effect of vibrational energy excited after many 
collisions (t^ n -t is large) in p’. Now, for example, assume that the 

rotational degrees are fully excited after two collisions or mean free 
paths and that the contributions of the vibrational energy to p' is 
zero 


tLi = 


P 


( 3 )( 2 )^= \ 


The estimates give meaning to the use made in the present paper of 
a p'/d of the order of unity for air undergoing slow and rapid flow 
changes. 
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APPENDIX B 


SOURCE FLOW WITH HEAT ADDITION 


Por "the general case of source flow with heat addition or variable 
total energy, the energy equation can no longer be used separately as 
was possible for source flow with constant total energy) rather, the 
equations of continuity, of motion, and of energy have to be integrated 
simultaneously. Furthermore, the amount of total energy addition or 
heat addition to the flow is no longer treated as a separable boundary 
condition, but all boundary conditions necessary to satisfy the system 
of differential equations have to be satisfied simultaneously. An 
exception is known to exist for the case of one-dimensional shock flow 
with heat addition (references 11, 12, and 19) where a separate integra- 
tion of the energy equation is made possible by use of certain values 
of the Prandtl number depending on the ratio M- ' /m- - This fact can also 
be seen from equation (7) in the present paper when F is taken constant. 

In the body of this paper the exponential deviations for terminated 
source flow from the main constant-total-energy solutions (curves 1 and 4) 
were discussed in some detail. Similarly, terminating deviations exist 
because of a heat source) however, as just indicated, for source flow 
the heat addition to the flow may no longer be treated as a separable 
boundary condition. In view of this complication, a much simpler problem 
is chosen to illustrate merely the basic nature of the effects caused by 
flow termination with a heat source. Assume that a wire is pulled with 
the speed u through a heat source at a fixed location. The differen- ^ 
tial equation for heat conduction "upstream" in the wire (the "downstream 
parts are heated as the wire is drawn through the source) is given by 


Integration indicates that a dying out of the heat-source effects occurs: 


where T M is the temperature of the unheated wire. 

In order to complete the physical picture, an estimate is made 
of the distance which the major effects, caused by the terminating 


dT _ _K_ ^fT 
dx _ pc p dx 2 


K d 3 ! 
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heat source, will penetrate "upstream"} the total effects will, of course, 
be conducted to 00 . In accordance with customary procedure (used, for 
example, for Prandtl's estimate of the shock thickness), the subtangent 
of the actual variation is taken as a measure of the penetration of the 
major effects. For the exponential temperature variation given in the 
preceding equation, the subtangent is 

K 

P uc p 

A rough estimate of this distance can be made by assuming that this 
measure holds approximately true for gas flow. If the various constants 
are left out from the pertinent expressions of the kinetic-gas theory 
(reference 6), the distance of penetration can be written as 

K ^ P ac p _ _z 
puCp ~ pucp M 

Attention is drawn to the fact that for this rough estimate in the 
neighborhood of Mach number of 1, a Mach number passed through by every 
shock, the penetration distance of the major effects of the terminating 
heat source is of the order of the mean free path. 
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APPENDIX C 

ENTROPY BALANCE 
General Development 


The discussion of the entropy balance is conveniently introduced by 
the complete energy balance for the same fluid element. The heat dQ 
added to the fluid element may be separated into the heat added from the 
outside by conduction dQ con ^ and the heat input due to internal genera- 
tion of heat by friction. The energy balance is expressed in terms of 
the familiar dissipation function 0 in the form (references 2, 20 , 
and 21) 


P 


dQ 

dt 


P 


dQ 

cond 

dt 


= 0 


(Cl) 


The term 


dQ, 


cond 


dt 


is known from the lav of heat conduction. The inflow 


or outflow of heat on one side of the element is 


-k grad T 

and the rate of heat addition by longitudinal heat conduction through 
two sides of the element (for source flow heat is not conducted trans- 
versally to the flow) is 


P 


^cond 

dt 


div(k grad T) 


For the purpose of subsequent comparison with the entropy balance, 
the energy balance is divided by T 


p " I div ( k § rad T ) = |r (C2) 

The entropy balance is made for an isolated system (references 22 
and 23). The entropy s of the fluid element is that of an open system 
since heat exchange by conduction occurs with the neighboring elements. 
The system can be made an isolated one by arranging at the boundaries 


Y 


NACA TN 2630 


4l 


of the fluid element an artificial entropy storage which compensates 
for the outflow of heat. The entropy flow through one wall of the fluid 
element is given by 


-k grad T 
T 

in contrast to the corresponding heat flow by conduction -k grad T. The 
entropy stores furnish for the entire fluid element the entropy increase 
rate 


div 


k grad T 
T 


Since, now the system is isolated, the second law may be expressed by 
the inequality 


P 


ds 

dt 


div 


k grad T ' > 


< 0 


(C3) 


To obtain the entropy balance in a more convenient form, the second 
term is developed with the aid of the general vector relation 


div(ab) = grad a.b + a div b 


where the vector b = k grad T and the scalar a 
entropy balance becomes 


1 

T' 


Therefore, the 


^ ds div(k grad T) k / , mX p \ 

p ^ T ) = 0 

or with the use of equation (C2) 

» ft - - I * T)2 6 0 


(C4) 


(C5) 


Since jZi ^ 0 (references 2, 20, and 21) and naturally k and T are 
positive, the correctness of the inequality is evident. 
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For the case that the total velocity of a steady flow is in the 

d s d s 

symmetry axis of the flow, the term p -jh becomes equal to pu — . 

d s 

Another form of the entropy balance may be obtained by expressing pu — 

dx 

in equation (C5) as a divergence which can be done by showing that 


pu = div(pvs) 
dx 


Specifically, 


div(pvs) = pv • grad s + s div(pv) 

The steady-flow continuity condition, however, expresses 

div( pv) = 0 


and, thus, 


div(pvs) = pv • grad s = pu 4^- 

dx 

Equation (C5) may be rewritten as 

aiv(p?s - * T ) > 0 

For source flow or for flow through slightly varying cross section 


or, since 


1 d „/ k dT\ - 

F dx l T dx/ “ 


puF = Constant, 


s 2 


JL_ > 
PUT dx/ 2 ~ 1 



(c6) 


where 0, 1, 2 represent subsequent points in the downstream 

direction. 
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For the case of shock flow with constant total energy in an infinite 
tube with constant cross section, it is not necessary to resort to this 
more complicated entropy balance. The reason is that, if in front (-«>) 
and behind (+ 00 ) the shock uniform conditions exist (that is, the velocity 
gradient and, thus, the temperature gradient are zero) , no heat conduc- 
tion will occur there, and thus, the shock as a whole may be regarded 
as an isolated system. This analysis is in essential agreement with a 
brief general statement in reference 19 concerning the reason why a 
negative entropy gradient occurring inside the shock (where a fluid 
element is an open system) does not violate the second law. The function 
of heat conduction in the entropy balance is, however, not directly 
mentioned in reference 19- 


It is well known that a certain isoenergetic mass flow with a given 
entropy can pass only through a certain minimum cross section. The 
physical reason for the penetration of an isoenergetic solution to cross 
sections smaller than the isentropic minimum of a sink flow, due to 
expansion in the continued sink flow, or in a curved minimum passage 
joined to the sink flow may be disclosed through the entropy balance, 
which is indirectly contained in the Navier-Stokes equation. The entropy 
balance is expressed by the inequality (C 6), where the subscript 0 would 
correspond to the isentropic minimum section and subscripts 1 and 2 
would refer to cross sections farther downstream. Since the problem is 
one of constant total energy, the entropies in the smaller downstream 
cross sections have to be progressively smaller than in the isentropic 
minimum section (references 2 and 2k). Thus, in order that the inequal - 


and, furthermore, both expressions have to be negative. Since for the 

present case of constant -energy flow ^ is the negative of — multi- 

dx 0 dx 

plied by a function of u, ^ has to be positive, and, in addition, it 

has to be sufficiently large as a consequence of the second law of 
thermodynamics. Since' the entropy balance is indirectly contained in 
the present solutions of the Navier-Stokes equation, the expansion 
curve 4 which penetrates beyond the isentropic minimum section has to 
meet these requirements of the entropy balance, provided, of course, 
that the Navier-Stokes equation is not used beyond the limit of its 
applicability. The entropy balance also indicates that a supersonic 
compression flow due to its positive temperature gradient will not be 
able to pass through a smaller than isentropic minimum section. 


Flow through Cross Sections Smaller than the 


Isentropic Minimum Section 


ity (C6) be satisfied, 



has to be greater than 
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APPENDIX D 

REDUCTION OF SOURCE-FLOW EQUATIONS TO SHOCK FLOW 
IN CONSTANT CROSS SECTION 


The case of a shock in constant cross section will be investigated 
with the aid of equation (l4) which is rewritten 


— - — f' 

, — P 5 

1 - 


21-tlu 2 - 1 

7-1 


byC^T 

7-1 


u 2 ^. - ^ 


7 - 1 
v 4tC p 


Vi - u 2 


U 


U I' 3 


(Dl) 


where 


1 - d log x 1 dx 
£ du x du 


Since for flow through constant cross section x is infinite, it is 

C 1 

eliminated from equation (Dl) by dividing by C 2 = which gives 

pux 


il = 

C 2 


il 


2 — — y u^ - 1 , 2 

7-1 I ' 


1 - u 2 C 2 — itZ_ u 2 /l - u2 C 2 2 

7-1 


( 7 - 1 

47 


/l - u 2 


u 


- C 2 u j^2 (D2) 


For the present flow x = oo and pu = Constant; thus 


ill = JJdxV 1 0 = 0 
C 2 2 x3\du/ 

Introducing the new variable 


pux. 


_|'_pudx_ r dx 
3 “ C 2 ~ c i du _ ^ du 
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gives 


V = 


u 


2-±4u2 

7-1 


1 - u 2 


JtZ — u 2 /l - u 2 
7-1 


(D3) 


The integration of equation (D3) is accomplished by transforming it in 
the following manner: Since 


d^L - u 2 = 


u du 
\J± - u 2 


du can be expressed as 


du = - j di/i 

u v 


—2 

u 


Now, equation (D3) can be written 



or 



d*n = 


71 d f 




]/l - u 2 


—2 27-1 

u + T| L 


47 jj3 


d \J± 


H 2 


7-1 


Multiplying the preceding equation by 


A - u 2 /i 


T] gives 


Vl - u 2 dr^ - d/l - u 2 T) 


y + 1 vT 2 - 


7-1 


u - 1 


u3 




u 2 d^l - u 2 


7 - 1 
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or, since 


a(yf 



u du 


-d 



7 + 1-2 




47 

7 - 1 



du 


From the integration of the preceding equation, the following equation 
re suits 



du + Constant 



7 - 1 

47 



du + Constant 


or 


£ u = J . + -._ u + ^1 ~ — = + Constant 
q 47 47 u 


or 


£ - u 2 


1, t 4 u + 2' ,. ~ 1 i + Constant 
47 47 u 
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Since in the final expression for constant -energy flow through the 

shock in constant cross section the velocity gradient — = — — 

J t\ C4 dx 

has to become zero, two infinities are to be expected for r\ as the 
uniform velocities u^ and u 2 in front of and behind the shock, 

respectively, occur. This condition requires an expression of the type 
(u - up)(u - u 2 ) or a quadratic form in the denominator. The quadratic 

form is obtained by multiplying and dividing the preceding equation by 
u which gives 


/l - u 2 u 

Z + - 1 u 2 + Z ~ 1 + Constant (u) 
47 47 


The denominator may also he written in the form 


X [^ 2 + y ~ + Constant (u)J 

Furthermore, across the shock, according to Prandtl, the relation 

O 

u-^ug = a* 


exists where a* is the critical velocity. Since u^ and Uq represent 
the ratios u 1 /u max and u 2 /u max 

7-1 

^ = 7 — 

and the denominator may be written as 


7 + 1 
by 



+ u l u 2 


+ Constant 
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Since 


(u - UjJ(u - ug) = u 2 - u (u-^ + u 2 ) + UjUg 

the expression inside the brackets represents (u - u-^u - u 2 ) and 
thus satisfies the conditions ahead of and behind the shock, if 
Constant = + u 2 ). The following expression, thus, may be written 

for x since T} = C^ 2^: 


_ 1 _ 47 

~ C k 7 + 1 


^1 - u 2 u du 
(u - U X ) (u - u 2 ) 


+ Constant 


This expression essentially agrees with the known expression for the 
velocity variation across a one -dimensional shock. (in comparing the 
result with that of reference 19 j note that ^' = 0 in the reference 
and that the variables used are not the same as those of the present 
paper . ) 
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TABLE I 


RATIOS OF 


ax 


GIVEN BY PLOTTED SOLUTION OF NAVIER -STOKES EQUATION TO BREAKDOWN 


Plot 


VALUES OF /— J 

\ ^ dx/Breakdown 


OF NAVIER -STOKES EQUATION 


[Plotted solutions are given in figures 1 to 3 3 


u 

rf , _ d log x 

1 / d log u\ 

/d log u\ 

/ p \ /d log u\ 

2^ du Id log x 1 

\ dx /Plot - breakdown 

d log u 

0 1 ” \d log xJ 

X 7 Plot 

Id log x) 
x Breakdown 

/ P \ /d log u\ 

2H \d iog xJ 

\ dx/Breakdown rioz 

Source flow) C 2 = 0.1) curve 1 

0.2 

0.04 

25.0 

17.50 

0.7 

.25 

.13 

7.69 

11.07 

1.44 

.3 

.27 

3.70 

7-57 

2.04 

A 

.90 

1.11 

4.09 

3.68 

• 5 

2.36 

.424 

2.47 

5.84 

.6 

7.17 

.140 

1.59 

11.38 

.65 

15.5 

.0645 

1.28 

19-9 

Source flow) C 2 = 0.1) curve 2 

0.1 

-0.92 

-1.09 

71.07 

-65.39 

.2 

-.33 

-3.03 

17.50 

-5A7 

• 3 

-.25 

-4.00 

7.57 

-1.89 

• 35 

-.30 

-3.33 

5.46 

-1.64 

A 

-.38 

-2.63 

4.09 

-1.55 

• 5 

-.70 

-1A3 

2A7 

-1.73 

.6 

-1.90 

-.526 

1.59 

-3.02 

Source flow) C 2 = 0.1) curve 3 

0.65 

14.1 

0.071 

1.28 

18.1 

• 7 

8.70 

.115 

1.04 

9.06 

•75 

2.38 

A20 

.84o 

2.00 

.8 

1.08 

.926 

.670 

.723 

.9 

.45 

2.22 

.384 

.173 



Sink 

flow) c 2 = -o.l) curve 4 

0.2 

-0.85 

-1.18 

-17.50 

14.9 

• 3 

-•70 

-1.43 

-7.57 

5.3 

A 

-.60 

-1.67 

-4.09 

2.46 

• 5 

-A9 

-2.04 

-2.47 

1.21 

.6 

-Ao 

-2.50 

-1.59 

.635 

• 7 

-.33 

-3.03 

-1.04 

.344 

.8 

-.26 

-3.85 

-.670 

.174 
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TABLE I.- Concluded 


RATIOS OF 




GIVEN BY PLOTTED SOLUTIONS OF NAVIER-STOKES EQUATION TO BREAKDOWN 


dx /Plot 

VALUES OF I— 


2 \i 


OF NAVIER-STOKES EQUATION - Concluded 


Breakdown 


[Plotted solutions are given in figures 1 to 3D 


u 

d log x 

1 _ U log u \ 

V \<J log X /Plot 

U log u\ 

/ P \ l \ 

ai du r log u ) 

\ dx /Plot __ \ los x ,fcreakdown 

\p ' — 

d log u 

\d lo S x / Breakdown 

/ P \ log u\ 

1^, du) \d log x/ plot 

\ dx / Breakdown 

Sink flow; Cg = -0.1; curve 5 

0.15 

1.40 

0.714 

-31.39 

-43.9 

.2 

.20 

5-00 

-17.50 

-3.50 

.3 

.18 

5.56 

-7.57 

-1.36 

.4 

.25 

4.00 

-4.09 

-1.02 

.5 

.35 

2.86 

-2.47 

-.866 

.6 

.50 

2.00 

-1.59 

-.794 

.7 

.76 

1.32 

-l.o4 

-.791 

.8 

1.18 

.847 

-.670 

-.790 

.9 

2.33 

.429 

-.384 

-.896 

• 95 

5.50 

.182 

-.247 

-1.36 

Sink flow; Cg = -0.1; curve 6 

0.1 

-0.02 

-50.0 

-71.1 

1.42 

.11 

-.06 

-16.7 

-58.7 

3.52 

.12 

-.n 

-9.09 

-49.2 

5.42 

.13 

-.17 

-5.88 

-41.9 

7.12 

.15 

-.68 

-1.47 

-31.4 

21.3 

.17 

-.89 

-1.12 

-24.4 

21.7 



Sink flow y C 2 = -0.1; curve 7 

0.46 

3.3, -3.13 

0.303, -0.319 

-3-00 

-9.89,9.38 

.5 

1.26,-1.05 

.794, -.952 

-2.47 

-3.12,2.60 


1.08, -.68 

.926,-1.47 

-1.97 

-2.13,1.34 

.6 

1.07, -.53 

•935,-1-89 

-1.59 

-1.70, .841 

.65 

1.13, -.v* 
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Curve 1 Expansion for infinite source flow 
Curve 2 Compression shock for infinite source flow 
Curve 3 Expansion due to change of boundaries 
" of source flow “ 



Figure 1.- Direction field and integration curves in plane of (f } plotted 
against u (u in log scale) for source flow with C2 = 0.1. (A large- 
size print of this figure is enclosed in the inside pocket of the hack 
cover page.) 
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Figure 2.- Direction field and integration curves in plane of plotted 

against u (u in log scale) for sink flow with C 2 = -0.1. (A large- 
size print of this figure is enclosed in the inside pocket of the hack 
cover page. ) 
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_ Curve 4 Expansion for infinite sink flow 
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Figure 3 -- Flow solutions in plane of u plotted against x (in log 
scale) for source flow with C 2 =0.1 (curves 1, 2, and 3) and sink 
flow with C 2 = -0.1 (curves 4 , 5 > 6^ and 7 )* (A large-size print of 
this figure is enclosed in the inside pocket of the hack cover page.) 
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(a) Direction field and integration curves in plane of 0' plotted 

against u (u in log scale) . 



(b) Flow solutions in plane of u plotted against x (in log scale). 

Figure 4.- Flow presentation when viscous effects of source or sink flow 

are zero, C2 = 0. 





Figure 5*- Detailed structure of direction field and integration curves 
in plane of 0' plotted against u (u in log scale) for C 2 = +0, 
the immediate neighborhood of C 2 = 0 for source flow. 
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Figure 6.- Detailed structure of direction field and integration curves 
in plane of 0' plotted against u (u in log scale) for C2 = -0, 
the immediate neighborhood of Cg = 0 for sink flow. 
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Curve 1 
Curve 2 
Curve 3 


Expansion for infinite source flow 
Compression shock for infinite source flow 
Expansion due to change of boundaries of source flow 



(a) Summary account of direction field and integration curves in plane 
of {' plotted against u (u in log scale) . 



(b) Flow solutions in plane of u plotted against x (in log scale). 

Figure 7*- Flow presentation of ~ the immediate neighborhood 

of C2 = 0 for source flow. 
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Curve 4 Expansion for infinite sink flow 

Curve 5 Compression shock for infinite sink flow 

Curve 6 Expansion due to change of boundaries of sink flow 




(b) Flow solutions in plane of u plotted against x (in log scale). 

Figure 8.- Flow presentation for C 2 = -0, the immediate neighborhood 

of C 2 = 0 for sink flow. 
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